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ABSTRACT 

We have identified four circumstellar disks with a deficit of dust emission from 
their inner 15-50 AU. All four stars have F-G spectral type, and were uncovered 
as part of the Spitzer Space Telescope "Cores to Disks" Legacy Program Infrared 
Spectrograph (IRS) first look survey of ~100 pre-main sequence stars. Modeling 
of the spectral energy distributions indicates a reduction in dust density by factors 
of 100-1000 from disk radii between ~0.4 and 15-50 AU, but with massive gas-rich 
disks at larger radii. This large contrast between the inner and outer disk has led 
us to use the term 'cold disks' to distinguish these unusual systems. However, 
hot dust [0.02-0.2 M moon ] is still present close to the central star (R <0.8 AU). 
We introduce the 30/13 /im flux density ratio as a new diagnostic for identifying 
cold disks. The mechanisms for dust clearing over such large gaps are discussed. 
Though rare, cold disks are likely in transition from an optically thick to an 
optically thin state, and so offer excellent laboratories for the study of planet 
formation. 

Subject headings: stars: pre-main-sequence — (stars:) planetary systems: pro- 
toplanetary disks 
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Introduction 



The evolutionary processes transforming massive, gas-rich circumstellar disks into tenu- 
ous, gas-poor debris disks are still not well understood. During this crucial interval, planets 
form and the remaining disk material is accreted or dispersed. The path by which this tran- 
sition proceeds remains uncertain. Do disks clear uniformly throughout via grain growth and 



settli ng; or is the process accelerated in the inner regions forming gaps (e.g. lAlexander et al. 



20061 )? Evidence of dust clearing should be visible in the infrared (IR) spectral energy distri- 
bution (SED). The Spitzer Space Telescope, with its wide wavelength coverage and increased 
sensitivity, is starting to reveal a new population of disks with unusual SEDs. 

Lack of mid-IR excess emission from disks has been interpreted a s a sign of dust clearing 
since the first disk SEDs were o bserved with IRAS (IStrom et al.lll989l ). Further examples and 
some spectra (e.g. HD 100546. iBouwman et al.ll2003l ) were obtained of a few Herbig Ae/Be 
stars with the ISO satellite. In this paper, we designate disks with this characteristic missing 
mid-IR emission as 'cold disks' due to the lack of emission from warm dust. These SEDs 
morphologically fall between the SEDs of classical optically thick disks, which show excess 
emission throughout the IR, and debris disks characterized by very weak far-IR excesses, 
leading them to be considered as transitional objects between the two classes. 

Spitzer spectra are particularly important to characterize the sharp rise, since they 
cover the critical 8-30 /im range. The earliest example discovered with Spitzer was CoKu 



Tau 4 modelled to have an inner hole within 10 AU (IForrest et al.l 12004!. iD'Alessio et al. 



2005). Three additional T Tauri stars, TW Hya, GM Aur and DM Tau flCalvet et al 



2002 



20051 ). have been identified as also containing cleared inner regions. One of the most exciting 
proposed explanations for these cold disks is that a planet has cleared a g ap in the disk 



and t hus that these inner holes may trace the presence of planetary systems (jVarniere et al. 
2006h . 



In this paper we present Spitzer spectra of four young stars of F-G spectral type that 
show a deficit of dust emission from the inner 15-50 AU of the disk yet strong excesses 
longward of 20 /im. These cold disk sources were identified in the first look survey of the 
Cores to Disks (c2d) Spitzer Legacy Project (jEvans et al.l 120031 ). comprising ~100 young (< 
10 Myr) stars, mainly T Tauri K and M stars, with circumstellar disks. All four sources 
are characterized by a steep, ~10x rise in flux beginning at 15 /im, indicating a sudden 
change in disk properties at a specific radius in a manner similar to that seen in the T 
Tauri transitional disks. Such objects are rare but provide an important window into disk 
evolution and planet formation processes. 
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Observations 



The 5-35 /im spectra were taken with the Infrared Spectrometer (IRS; iHouck et al. 



2004f ) on Spitzer using both low (R=160, As L :5.2-14.5/im, A LL :14.0-38.0/im) and high (R=600 
As#:9.9-19.6/im, A^: 18.7-37 .2/xm) resolution modules. Spectra were extracted from the 
Spitzer Science Center (SSC) Basic Calibrated Data (BCD) images, generated by pipeline 
S13. For the low resolution spectra, the SSC pipeline full aperture extraction was used. 
For the high resolution modules, the c2d extraction, based on a combined sine fitting of the 



spect ral trace to account for bad pixels and background emission, was used (see lLahuis et al. 
20061 for further detailsQ). 



MIPS SED spectra, taken as part of the c2d project, are included in the SEDs for all 
four cold disk sources. For each source the BCDs were coadded using MOPEX. The coadded 
images were extracted with IRAF, using an optimized extraction with a 3 or 5 column wide 
aperture. 



3. Distinguishing characteristics 

The large (~100 object) c2d IRS first look sample allows comparison between the four 
cold disk sources and the remaining systems to identify trends which might be diagnostic of 
their evolutionary state. The small number of cold disk sources (< 5% identified), both in 
the c2d sample and in the literature, suggests that this condition is rare either due to rapid 
evolutionary timescales undergone by all stars or due to an unusual condition unique to a 
small sample of stars. 

The four cold disk sources are clearly differentiated from the majority of the c2d 
star+disk systems by the ~10 fold increase in IR flux between 10 and 30 /an. To char- 
acterize this increase in dust emission, we have used the 30/13 /im flux ratio (see Figured]). 
These wavelengths have been chosen to avoid strong silicate features, but include the full 
increase in continuum emission. The majority of the sample has emission that increases by 
a factor of 2.3 ± 1.4 between 13 and 30 /im while the transitional disks rise by factors of 
5-15. 

Interestingly, out of a sample comprised of predominantly low-mass stars, the cold disk 
sources are all of intermediate mass with spectral types of F and G. The majority of K and M 
stars and the higher mass A and B stars have low 30/13 /im ratios. However, there are only 



1 Lahuis et al. 2006 is available at: 



http://ssc.spitzer.caltech.edu/legacy/c2dliistory.html 
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a handful of A and B stars in our sample so it is difficult to draw any significant conclusions 
about such stars. 

All four cold disk sources show polycyclic aromatic hydrocarbon (PAH) features, partic- 
ularly at 11.3 A tm. Such emission is uncommon in the c2d sample with only ~10% displaying 
PAH features ( iGeers et al.l 12006k see their Figure 5 for blow-ups of the PAH bands). The 
presence of PAH features may be enhanced in these cold disks because the lack of dust in 
the inner disk lowers the mid-IR continuum flux creating a stronger line-to-continuum ratio, 
thus facilitating detection. 

The four disks are also characterized by weak to non-existent 10 fim amorphous silicate 
features. LkHa 330 is the only one that shows an unambiguous, though low contrast, 10 
/mi feature. The spectrum of HD 135344 includes the wavelength region of the full 10 fim 
band but shows no silicate feature. The spectra of SR 21 and T Cha begin at 10 /iin, but 
appear to have only weak silicate emission, if any. The broader 20 //m silicate feature is 
harder to isolate from continuum dust emission, particularly with the sharp rise in the SED 
beyond 15 fim. However, some 20 fim emission does seem to be present. This feature traces 
colder regions of the disk than does the 10 fim feature, and thus indicates the presence of 
amorphous silicates further out in the disks. 

The Ha equivalent width is often used as a tracer of accretion. Different dividing lines 
between non-accreting weak-line T Tauri stars (wTTs) and classical T Tauri stars (cTTs) 
have been proposed bu t it makes little dif ference to the classifi cations in these cases. SR 21 



(0.54 A in absorp tion, Martin et al.lll998h and T Cha (2-10 A, iGreeorio-Hetem et alJll992 



Alcala et al.lll995l ) are nomin ally wTTs, although T Cha is highly variabl e and close to the 
cutoff. LkHa 330 (11-20 A, iFernandez eiTall Il995l . ICohen fe Kuhl Il979h and HD 135344 
(17.4 A. lAcke et aDl2005h are clearly cTTs. 



4. Modeling 

Modeling is necessary in order to interpret the SEDs in terms of the physical structure 
of the disk. In particular, can a disk model with a gap accurately reproduce the SEDs, 
including the steep rise between 13 and 30 /mi? Is such a fit possible without resorting to a 
gap? 



The disks were modeled with the 2-D radiative transfer code RADMC ( iDullemond &: Dominik 



20041 ). The steep factor of 10 rise in flux between 13 and 30 /mi prevents these disks from 
being fit well by conventional disk models. We introduce a very wide gap with an inner 
radius, -RGap,in, near 1 AU and an outer radius, -Rcap.out, near 30 AU (see Tabled] for specific 
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values). In order to model the steep change in emission associated with i?Gap,out, grid refine- 
ment is introduced at -Rcap.out as well as the inner dust rim, i?Disk,in- A physical reduction 
in dust density is only one possible scenario which could result in this SED shape (see §5). 

Input parameters for these models include the stellar mass, M*, radius, R*, and effective 
temperature, T e g. Kurucz models are used for the stellar photospheres. Where p ossible, val- 



ues fo r stellar and disk properties are taken from the literature (see Table [p. The lSiess et al. 



(120001 ) pre-main sequence stellar tracks were used to check that the parameters were consis- 
tent. The effects of modest differences in the stellar properties on the mid-IR portions of the 
spectra are small. For example, a 100 K change in the stellar temperature results in a change 
in -RGap,out of 2-3 AU. Optical and (sub)m illimete r photo metry is often not simultaneous and 
is dereddened using the extinction law of iDraind (120031 ). 



The disks are assumed to be flared with surface height H such that H/R oc i? 2//7 as in 
Chiang fc Goldreichl (119971 ). The models do not calculate the hydrostatic equilibrium self- 
consistently and here the pressure scale height is set at the outer disk edge, -Rdisk, which 
is assumed to be 300 AU for all the models. H/i? disk values slightly lower than predicted 
by hydrostatic equilibrium are needed to fit the outer disks of both HD 135344 and T Cha, 
which produces a smaller silicate feature even in the models with no gap. This could indicate 
dust settling in the out er disk, which mig ht be expected from the older ages (8-10 Myr vs 
2-3 Myr) infered by the ISiess et al.l (120001 ) tracks. The dust composition is set in all models 
to a silicate: carbon ratio of 4:1. The grain sizes range from 0.01 to 10 fim with a power-law 
index of -3.5 and a du st mass of Mrw^ma ii- The upper size limit of 10 /xm is large compared 
to interstellar grains (IMathis et al.lll977l ). but was necessary to account for the lack of any 
10 (Urn amorphous silicate feature. This indicates a skewing of the size distribution likely 
through grain growth. No crystalline silicates are included. A midplane layer of 2 mm grains 
containing 90% of the total dust mass, with the rem aining 10% in Mnnst^ man, has been 



added to account for the (sub) millimeter slope (as in iDullemond fc Dominiki 120041 ). The 
inner edge of the disk, -RDisk,in, was set at the radius where the dust sublimation temperature 
of Tdust ~ 1500 K is reached. A slightly puffed up inner rim was included with if/ i?Disk,in 
being set at 0.03, just above the hydrostatic equilibrium value. The disk inclinations are not 
well constrained so i ~ 30° was used arbitrarily f or SR 21 and LkHa 330. T Cha is more 
edge on so i ~ 75° wa s adopted (lAlcala et al.lll993l ) and HD 135344 is more face-on with i ~ 



10° (iDent et al.ll2005l ). However, until the disk is nearly edge on and begins to obscure the 
central star, varying the inclination produces little effect on the gap properties. 

The optimal size of the gap necessary to fit the mid-IR is found by minimizing the \ 2 
of the fit (see Figure [3]). -Rcap.out is tightly constrained by the steep rise at 13 /zm, while 
-RGap,in is less well constrained. The dust density in the gap was reduced by 10~ 6 compared 
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to the equivalent models with no gaps, and the minimum reductions needed were factors of 
10-100. 

All of our sources have 1-10 /im excesses which require 0.02-0.24 M moon of dust between 
-Roisk.in and -Rcap.in, a fraction of only 10~ 7 to 10~ 6 of the total refractory dust mass available. 
This matter is optically thick to the stellar light at the midplane. If we remove all matter 
from the disk inward of -Rcap.out, then the near-IR fluxes are underpredicted. These excesses 
clearly show that at least some dust must be close to the star. We model the 1-10 /im excess 
by keeping the inner disk intact and only inserting a gap from -RGap,in to -RGap,out, where 
-RGap,in > -Roisk.in- For T Cha, this results in a particularly small, hot dust region inside the 
gap which would have a very short lifetime, of order thousands of years at most, without a 
continual influx of material. 



Discussion 



We have identified four cold disks around F and G type stars with unusually steep flux 
rises between 10 and 30 /zm, whose SEDs can only be modelled with wide gaps of inner radii 
of 0.2-0.8 AU and outer radii 15-50 AU. These gaps are generally larger than those inferred 
for the 4 T Tauri stars with similar 30/13 /im ratios found so far, which have outer gap 
radii of 10-24 AU. Another difference is that our sources have small 1-10 /im excesses which 
demand that hot dust exists between the inner edge of the gap and the star, i.e. our sources 
have gaps rather than holes. The statistics on sources of different spectral types are still too 
small to conclude whether this is a general trend or peculiar to our sources, and whether 
this extends to Herbig Ae/Be stars. 

Dust clearing related to planet formation would be one of the most exciting explana- 
tions for the observed SEDs. There is also the possibility of the disk being disrupted and 
the dust cleared by a stellar companion. HD 135344 and SR 21 are part of wide binarie s 
with separations of 20'/4 and 6'.'4, respectively (jCoulson fc Waltherlll995l ; iPrato et al.ll2003l ). 
Although no companions are currently known within the modeled -Rdisk; close binaries (<50 
AU) cannot be ruled out for any of the sources. 

An other proposed sc e nario for quickly clear ing the inner disk region is photoevapo- 
ration (IClarke et al.l l200ll . Alexander et al.l 120061 ) . This physical process occurs when the 
photoevaporation rate driven by the ionizing flux from the central star matches the viscous 
accretion rate, resulting in an inner hole. The predicted size of the inner hole is given by 
R g = GM^/c^, with T ~ 10 4 K to give c s ~ 10 km/s, leads to predicted hole radii of 13-18 
AU for 1.5-2 M stars, although a more rigorous examination of the gas dynamics revealed 
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that this equation over-estimates the inner hole size (jLiffmanl 120031 ) . However, models of 
Herbig Ae stars (M*=2.5 M ) that combine photoevaporation with viscous evolution and 
differential radial motions of dust and gas predict that t he inner disk clears quickly but leaves 
gas-poor dust rings at 10-100 AU (ITakeuchi et al.ll2005l ). The accretion rates for all four cold 
disks, with the possible exception of SR 21, are too high to make this scenario likely. 

An alternative to physically removing the dust is to let it grow beyond the size at 
whi ch it efficiently radi ates as a blackbody so that it no longer emits strongly in the mid- 
IR (ITanaka et al.l 120051 ) . T here is general evidence for grain growth in d isks from both 
mid-IR and millimeter data (IKessler-Silacci et al.ll2006l . iRodmann et al.ll2006l ). The lack of 
strong 10 /im amorphous silicate features also points to grains having grown beyond inter- 
stellar si zes. There has been m u ch rec ent modeling work on clearing disks through grain 
growth. iDullemond &: Dominiki (120051 ) found that cold disk SEDs could be produced by 
dust coagulation but the timescales were too fast. Replenishment of the dust was necessary 
to match observed lifetimes of the disks. If replenishment processes such as fragmenta- 
tion occur preferentially i n the innermost region due to higher temperatures and densities 
(iKenyon &: Bromley! 12004 ) , this migh t be su fficient to produce the fraction of a lunar mass 
needed close to the star. iRice et al.l (120061 ) invoke dust filtration by an embedded planet 
whereby large grains pile up at the outer edge of the gap while small grains and gas pass 
through, thus accounting for both the hot small dust grains needed to fit the near-IR and 
the high mass accretion rates. The high spatial resolution of ALMA should be able to test 
these scenarios by searching for and imaging the emission from (sub)millimeter to centimeter 
sized grains. 
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Table 1. Model parameters 



Source 


Spectral 


A v 


Distance 


M, 


T e ff 




-RDisk,in 




-*^Gap,out 


■Mrjust, small 




-ffp(-Rdisk)/ 




Type 


(mag) 


(PC) 


(M e ) 


(K) 


(L Q ) 


(AU) 


(AU) 


(AU) 


(10- 6 Mq) (M lunar ) 


-Rdisk 


LkHa 330 


G3 a 


1.8 Q 


250 b 


2.5 C 


5800 


16 


0.27 


0.8 


50 


5 


0.24 


0.15 


SR 21 


G2.5 d 


Qd 


160 


2.5 d 


5800 


24 


0.25 


0.45 


18 


15 


0.10 


0.17 


HD 135344 


F4 e 


0.5 f 


84 e 


1.8 


6600 e 


6.8 


0.18 


0.45 


45 


5 


0.10 


0.13 


T Cha 


G89 


1.5 s 


66^ 


1.5 9 


5600 


1.4 


0.08 


0.2 


15 


3 


0.025 


0.11 



1 a - ICoheri fc Ktihil Jl979h. b - lEnoch et all d2006h. c - lOsterloh & Beckwithl Jl995h. d - 
Dunkin et all |l997t ). f - iMalfait et al.l |l998t )~q - lAlcala et al.1 |l993l ). h - Iwichmann et al. 
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Fig. 1. — Spitzer IRS 30/13 /xm flux ratios for the c2d first look disk sample, along with 
those of DM Tau, GM Aur and CoKu Tau 4. The cold disks have much larger 30/13 /xm 
ratios than does the majority of the sample. Outliers are labeled. 
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Fig. 2. — Spitzer IRS and MIPS spectra of the four cold disk sources reported here. The 
dashed line is the best fit star+disk model with a gap (see Table 1) and the dotted line is the 
equivalent model with no gap. JHK photometry is from 2MASS. IRAS photometry and our 
Spitzer photometry is shown in the mid- to far-IR. Optical and submillimeter photometry 



is from the literatur e ( Alcala et al 



1993 



Henning et all 1 19931 : lOsterloh fc Beckwithlll995h . 



Andre &: Montmerk 1994: Fernandez et al. 1995 
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Fig. 3. — Chi squared contour maps of 2-D disk model fits to the SEDs with different inner 
and outer gap radii (clockwise from top left: LkHa 330, SR 21, T Cha and HD 135344). 



